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We summarize the main characteristics and recent results on _B Xu(.vi decays of a model based on soft-gluon 
resummation and an analytic time-like QCD coupling. 



1. Introduction 

By comparing various spectra in the semilep- 
tonic inclusive B decays B I , where t is 

a fixed lepton species = e, //) and Xu the frag- 
mented u quark, with the predictions of a model 
including non-perturbative corrections to soft- 
gluon dynamics through an effective QCD cou- 
pling [T], we obtain a value for the \Vub \ Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element |2l3j 
\Vuh\ = (3.76 ± 0.13 ± 0.22) x IQ-^ being the 
errors experimental and theoretical, respectively. 
The model involves the insertion, inside standard 
threshold resummation formulae, of an effective 
QCD coupling a5(fc^), based on an analyticity 
requirement and resumming absorptive effects in 
gluon cascades [4]. By construction, Q!s(/e^) has 
no Landau pole and saturates at small scales: 

hm as{e) = 1 « 0{l). 

This model, which has no free parameters, de- 
scribes B-meson fragmentation data at the Z'^ 
peak rather well [5] , where — unlike B decays — 
accurate data are available and there is no uncer- 
tainty coming from the CKM matrix elements. 

In the following, we describe the phenomelogi- 
cal model used, then the extraction of \Vub\, and 
finally the results, followed by the conclusions. 



2. Threshold resummation with an effec- 
tive coupling 

Factorization and resummation of threshold 
logarithms in semileptonic decays leads to an ex- 
pression for the triple-differential distribution, 
the most general distribution, of the following 
form [S]: 

1 d^r 

, , , = C[x,w;asiQ)]a[u-Q] 
i axawau 

+D[x,u,w;as{Q)] , 

where: 

2Ee Q 1 - - (2mx/Q)2 

X iv — U — — 

"lb ' TTib" 1 + ^1 - (2mx/Q)2 

and the hard scale Q — 2Ex, with Ei, Ex and 
mx being the charged lepton energy, the total 
hadron energy and the hadron mass, respectively. 
F = r(as) is the inclusive width of decay B — > 
Xu £ ve- Furthermore, C[x,w;as] is a short- 
distance, process dependent hard factor; a[u; Q] is 
the universal QCD form factor for heavy-to-light 
transitions, resumming to any order in as the se- 
ries of logarithmically enhanced terms to some 
logarithmic accuracy; D[x,u,w;as] is a short- 
distance, process dependent, remainder function, 
vanishing in the threshold region u — > and in 
lowest-order in as- The heavy fiavor decay form 
factor has an exponential form in Mellin moments 
iV-space [7]. We apply a change of renormaliza- 
tion scheme for the coupling constant as as- 
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The QCD form factor (j[u; Q] has been numer- 
icahy computed for different values of as{mz) 
in[T]. 

Threshold suppression — the main theoretical 
ingredient for the measurement of | Vut \ — is rep- 
resented by the factor 



W{a,b) 



_T[B^Xuiiye,pGia,b)] 



T[B 
1 d^T 



a<ip{x ,w ,u) <b 



r dxdwdu 



dxdwdu < 1 . (1) 



where W{a, b) — 1 when integrated over all the 
phase-space of a;, w and u. For threshold re- 
summed spectra of i? ^ Xu £ vi decays at next- 
to-leading order see Ref. [8]. 

3. \Vub\ extraction 

Experimentally, given a kinematical variable p, 
such as for example the energy i?^ of the charged 
lepton, one measures the number of -B's decaying 
semileptonically to X^ with p in some interval 
(a, 6), divided by the total number of produced 
B's (decaying into any possible final state): 



B[p(^{a,b)] 



N\B 



Xuivi, p e (a, 6)] 



N[B ^ (anything)] 
This branching ratio can be written as: 

BsL 



B[pe{a,b)] 



1 



W{a,b), 



(2) 



(3) 



where we have defined the semileptonic branching 
ratio: 



BsL 



T{B ^ X,£iy,) + T{B ^ X^iue) 
T[B ^ (anything)] 



and the ratio of (6 c)/{b u) semileptonic 
widths: 



7^, 



c/u 



T{B Xc£iy{) 
TiB ^ X^li^i) 



(4) 



Since Bsl is rather well measured, we use the ex- 
perimental determination [9]: 

Bsl = 0.1066 ± 0.0020. 

With this method there is no ml dependence 
(with the related uncertainty) which would ap- 
pear using the theoretical expression of r(i? 



Xu (.I'i), because one has to compute only the ra- 
tio of widths TZc/u 
The semileptonic b 



tio of widths TZc/u and not the absolute widths. 



width is written as: 



T{B Xcive) 



G|m5 iKbl 



1927r3 
xI{p)F{as)G{as.p), 



(5) 



where p = ^ w 0.1. The function I{p) ac- 
counts for the suppression of phase-space because 
of mc ^ lOJ : 

I{p) = 1 - 8/9 + 12p2 log i + _ p4 _ 

P 

Note that there is an (accidental) strong depen- 
dence on the charm mass mc^ because of the ap- 
pearance of a large factor in the leading term in 
p, namely — 8. As far as inclusive quantities are 
concerned, the largest source of theoretical error 
comes indeed from the the uncertainty in p. Most 
of the dependence is actually on the difference 
mtj — rTic, which can be estimated quite reasonably 
with the Heavy Quark Effective Theory (HQET). 
Finally, the factor G{as,p) contains corrections 
suppressed by powers of as and p: 



G{as,p) 
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withG„(0) = 0. Note that G(0,p) = G(as,0) = 
1. By inserting the above expressions for the 
semileptonic rates, one obtains for the perturba- 
tive expansion of TZc/u- 

TZc/u = TZc/u{p,asA^ub\/\Vcb\) = 



\Vcbl 



/(p)G(a5,p). 



This method actually provides a measurement of 
the ratio |Ki6|/|V^cfc|, but since the error on \Vcb\ 
is rather small and theoretically well understood, 
one is basically measuring |Ki6|I3 

3.1. Quark masses 

The approach we are following to compute \ Vub\ 
can be subdivided into two parts: in the first 



^ The average of determinations of \Vcb\ coming from a 
global fit to the B — > Xcdui and b — > 57 moments in 
the kinetic and IS schemes, in good agreement with each 
other, is |V;6|= (4f.6±0.6) X 10"^ [TTE]- 



3 



part we compute the triple-differential distribu- 
tion and in the second part the \ Vub\ value. As far 
as the triple differential distribution is concerned, 
as the whole process is described in a perturba- 
tive framework, we do not distinguish between 
the mass of the B meson and the pole mass of the 
b quark, i.e. we consistently assume mb = tub- 

Once we compute the ratio |K,f)|/|V"c6|, we use 
the standard HQET formulas, based on the b~ 
quark mass. Thus, the 6-quark mass is intro- 
duced in our formulation of the ratio. Since 
quarks are confined inside observable hadrons, 
their masses cannot be directly measured and 
their values are biased by the selected theoretical 
framework. We have performed the calculation 
in the MS mass scheme. The MS masses for the 
b and the c quark are taken rnbirrib) = 4.243 ± 
0.042 GeV and rn^im^) = 1.25 ± 0.09 GeV [TTT9] . 
respectively. However, we have considered the 
pole-mass scheme as well in order to take into 
account the uncertainties coming from a different 
scheme definition. 

4. Results 

We calculate | Vub \ for all the experimental anal- 
yses. They are categorized according to the kinc- 
matical distribution looked at, where selection 
criteria are applied to define the limited phase- 
space on which the branching ratio is computed. 
They are: the lepton energy {E^), the invariant 
mass of the hadron final state (Mjf), the light- 
cone distribution {P^= Ex — \px\, Ex and px 
being the energy and the magnitude of the 3- 
momentum of the hadronic system), a two di- 
mensional distribution in the electron energy and 
s™"^, the maximal Mx"^ at fixed and E^. 
Moreover, as described in Ref. [1], we will look 
only at the range where data are not affected by 
potential b ^ ctvg background: 2.3 GeV< Eg< 
2.6 GeV, although the method works for lower 
lepton energy as well. 

We compute \Vub\ for each of the analyses 
starting from the corresponding partial branch- 
ing fractions. Then, we determine the average 
\Vub\ value using the HFAG methodology [TT]. 

Table [1] reports the extracted values of \Vub\ 
for all the uncorrelated analyses and their corre- 



sponding average. The errors are experimental 
(i.e. statistical and systematic) and theoretical, 
respectively. The average is: 

V^b = (3.76 ± 0.13 ± 0.22) x 10"^ 

consistent with the measured value of \Vub\ from 
exclusive decays |11| and the indirect measure- 
ment [T2l. 

The table shows also the criteria used for 
the determination of the partial branching ra- 
tio {AB). The theoretical errors are considered 
completely correlated among all the experimental 
analyses, when performing the average. The \Vub\ 
values and the corresponding average are plotted 
in Figure [TJ 



2 3 4 

\YJ [X 10-^] 

Figure 1. \Vub\ values for the uncorrelated anal- 
yses and their average. 

Several sources of theoretical errors have been 
considered and shown in Table [2l another 
method, based on the absolute values of the decay 
rates, to extract the value of |T4b|, the pole in- 
stead of the MS mass scheme, the order at which 
the rate is computed from the exact NLO to the 



BABAR (E,) 
3.46 ±0.14 ±0.23 
BELLE (E,) 

3.25 ±0.17 + 0.22 -0.21 ' ' ' ' ' 
CLEO (EJ 

3.49 ± 0.20 ± 0.23 " — ► 

BABAR 

4.04 ±0.19 ±0.24 

BELLE 

3.93 ± 0.26 ± 0.23 
BABAR m^-q^ 
4.14 ±0.26 ±0.23 
BELLE sira. ann. (iTijj, q") 
3.95 ± 0.42 ± 0.22 
BABAR (E„ s™") 
3.87 ± 0.26 + 0.23 - 0.24 
BABAR P+ 

3.45 ±0.22 + 0.21 -0.37 ' ' ^ 

Average +/- exp + theory - theory 
3.76 ±0.13 ±0.22 
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Table 1 

The first column in the table shows the uncorrelated analyses, the second column shows the corresponding 
values of \ Vub\, and finally the last column shows the criteria for which AB is available. The final row shows 
the average value of \Vub\- The errors on the \Vub\ values are experimental and theoretical, respectively. 
The experimental error includes both the statistical and systematic errors. 



Analysis 


iKbKio-^) 


AB criteria 


BaBar (Ee) [l3] 
Belle (Ee) [H] 
CLEO (Eg) [H] 
BaBar (Mx) [S] 
Belle (Mx) 
BaBar ((Mjf,g2)) [16] 

Belle iMx,q^)M 
BaBar (E^f, s""'^)[19] 
BaBar (P+) [IS] 


3.46± 0.14 to 'it 
3.20± 0.17 tin 
3.49± 0.20 tin 
4.04± 0.19 tlit 
3.93± 0.26 tofa 
4.14± 0.26 tin 
3.95± 0.42 tin 
3.87± 0.26 tin 
3.45± 0.22 tin 


Ei> 2.3GeV 
Ei> 2.3GeV 
Ei> 2.3GeV 
Mx< 1.55 GeV 
Mx< 1.7 GeV 
Mx< 1.7 GeV, >8 GeV^ 
Mx< 1.7 GeV, g2 > 8 GeV^ 
Ei,> 2.0 GeV, s™°^< 3.5 GeV^ 
P+< 0.66 GeV 


Average 


3.76± 0.13 tm 





approximate NNLO [5U] , the variation of the pa- 
rameters in the computation of \Vub\ within their 
errors, as given by the PDG [9]. What we cannot 
change is the modelling of the threshold region, 
which is fixed in our model, because it has no free 
parameters. That is the factor W{a,b). The er- 
ror on the modelling of the threshold region can 
only be estimated in an indirect way, by consider- 
ing different decay spectra, in which presumably 
threshold effects enter in different ways. Even 
though our model on soft-gluon dynamics is for- 
mally without free parameters, we may say that 
we have constructed it, out of many possibilities, 
as a kind of "expert system". Once "trained" 
by giving beauty fragmentation data in input, it 
should predict reasonable beauty decay spectra. 

Table O shows the | Vub \ averages for different 
analysis categories. Note that the {Mx,q^) anal- 
yses tend to have the largest values of \ Vub\, while 
the endpoint analyses the smallest. The larger 
value of \Vub \ coming from the analysis of the dou- 
ble distribution in {mx,q^) is expected on quali- 
tative basis. The lower cut on therefore signif- 
icantly reduces the hard scale Q from the "nat- 
ural" value Q = iTiB, where our model has not 
been tested as it has been constructed to describe 
B-decay spectra having the (maximal) hard scale 
Q = ms- 



Finally, using the value of AB corresponding 
to the lowest lepton energy cut for the endpoint 
analyses, the value of \Vub \ is 3% higher than what 
we quote adopting a cut at 2.3 GeV, qualitatively 
consistent with the expectactions due to the ob- 
servation of a larger number of events in that re- 
gion than predicted [T . 

5. Conclusions 

We have analyzed semileptonic B decay data 
in the framework of a model for QCD non- 
perturbative effects based on an effective time-like 
QCD coupling, free from Landau singularities. 

Our inclusive measurement of the \Vub\ CKM 
matrix element is: 

\Vub\ = (3.76 ± 0.13 ± 0.22) x 10"^ 

The errors on the \Vub\ value are experimental 
and theoretical, respectively. The experimental 
error includes both the statistical and system- 
atic errors. This value is fully consistent with 
the determination from exclusive decays (|V^6| = 
(3.51 ± 0.2ll° j^) X 10-3) El and from an indi- 
rect estimate (jKbl = (3.44 ±0.16 x IQ-^), 
whilst other methods show a discrepancy up to 

^The value is obtained from an average of the Lattice QCD 
determinations 
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Table 2 

The first column of the table shows the different contributions to the theoretical errors, the second column 
shows the corresponding variation, and finally the third column shows the percentage contribution with 
respect to the \Vub\ value. 



Theoretical Errors 


Contribution 


Variation 


Error (%) 


as 

\Vcb\ 

rub {GeV) 
mc (GeV) 

\Vub\ method 
pole mass (GeV) 
approx. NNLO rate 


0.1176 ±0.024 
(41.6 ±0.6) X 10-3 
4.20 ±0.07 
1.25 ±0.09 
0.1066 ± 0.0020 

4.7 <mb< 5.0, 1.47 < < 1-83 
3.34 < nib - rric < 3.41 


±0.6 ^ 3.5 
±1.4 
±0.6 
±4.4 
±1.0 
+0.8 

-1.3 ^ -5.2 

+2.0 



Table 3 

The table contains the \Vub\ values for several analyses and the corresponding averages. The errors on 
the \Vub\ values are experimental and theoretical, respectively. The experimental error includes both the 
statistical and systematic errors. 



\Vub\ for endpoint analyses (10 



BaBar (Ee) [13] 
Belle (El) [g 
CLEO (Ee) [H] 


3.46± 0.14 toM 
3.20± 0.17 toil 
3.49± 0.20 t|];23 


Ei> 2.3 GeV 
Ei> 2.3 GeV 
Ee> 2.3 GeV 


Average 


3.42± 0.15 toi^i 






Vub\ for Mx analyses (10 


BaBar (Mx) [IS 
Belle (Mx) 


4.04± 0.19 t'^aiil 
3.93± 0.26 tall 


Mx< 1.55 GeV 
Mx<1.7 GeV 


Average 


4.00± 0.16 toilt 




\Vub\ for {Mx,q') analyses (IQ-^) 


BaBar (Mx,g^)[T6] 
Belle {Mx,q^)m 
Belle {Mx,q^W, 


4.14± 0.26 toifs 
4.21± 0.37 toll 
3.95± 0.42 tall 


Mx< 1.7 GeV, g2 > 8 GeV^ 
Mx< 1.7 GeV, >8 GeV^ 
Mx< 1.7 GeV, g2 > 8 GeV^ 


Average 


4.13± 0.21 tall 






Vub\ for P+ analyses (10^^) 


BaBar (P+) [16] 
Belle (P+) [n] 


3.45± 0.22 toir 
3.73± 0.32 t^all 


P+ < 0.66 
P+ < 0.66 


Average 


3.55± 0.19 tals 





« +2 a J+j from those. We argue that the main 
difference of our model with respect to previous 
ones is a smaller suppression of the threshold re- 
gion. 
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